Adiabatic quenches through an extended quantum critical region 
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By gradually changing the degree of the anisotropy in a XXZ chain we study the defect formation 
in a quantum system that crosses an extended critical region. We discuss two qualitatively different 
cases of quenches, from the antiferromagnetic to the ferromagnetic phase and from the critical to the 
antiferromegnetic phase. By means of time-dependent DMRG simulations, we calculate the residual 
energy at the end of the quench as a characteristic quantity gauging the loss of adiabaticity. We 
find the dynamical scalings of the residual energy for both types of quenches, and compare them 
with the predictions of the Kibble-Zurek and Landau-Zener theories. 
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When, by changing some external parameters, a quan- 
tum system crosses a phase transition at a finite speed 
it is unable to reach its equilibrium (or ground) state no 
matter how slow the rate of the quench is. The reason 
is that the anomalous slow dynamics close to the criti- 
cal point prevents the system to follow adiabatically the 
external drive. As a result, supposing that the dynamics 
takes the system to a symmetry broken phase, a num- 
ber of defects will appear in the final state. Dynamical 
defect formation has important implications in a wide 
spectrum of problems ranging from the study of phase 
transitions in the early universe P, Q and in superfluid 
systems 0, 0, Hi , to quantum annealing @, 0, d, 0, 
and adiabatic quantum computation [ll|. 

Early works on the adiabatic crossing of a phase transi- 
tions dealt with classical system where the external con- 
trol parameter is the temperature. This problem was 
more recently explored also in the case of a quantum 
phase transition l3, 13] 
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tense theoretical activity ( 

m, m [H, [H, m, [H, and references therein). The 



large body of results obtained so far for the adiabatic 
crossing through a critical point are in agreement with 
the Kibble-Zurek (KZ) theory. At the roots of the KZ 
mechanism there is the hypothesis that the dynamics of a 
system close to a continuous phase transition can be con- 
sidered adiabatic or impulsive depending on the vicinity 
to the critical point. The determination of the point in 
which the system stops following the external drive leads 
to the determination of the scaling of density of defects 
and other observables as a function of the quench rate. 
In the limit of very slow quenches, defect formation can 
be also understood by means of the Landau-Zener (LZ) 
theory applied to the ground and the first excited states. 

The passage through a critical point is not the only sit- 
uation one can envisage in this context. Another paradig- 
matic case is when the evolving many-body system is 
quantum critical in an extended region of the parameter 
space. This is the topic of the present work. The sys- 
tem we use to illustrate this situation is described by the 



cross the different phases is the anisotropy coupling. The 
study of adiabatic quenches in the XXZ model allows to 
test several aspects of the problem which could not be ad- 
dressed previously. The system has an extended critical 
line instead of a critical point (a similar issue was recently 
considered in for the Kitaev model). Moreover, the 
boundaries of the critical region are characterized by dif- 
ferent exponents, hence one can test if the defect density 
is controlled by what happens before or after the pas- 
sage through the critical point [2§|. Finally, this model 
allows to study defect formation in the presence of dy- 
namical constraints. Because of the conservation of the 
total magnetization, it will be impossible for the system 
to reach the local ground state no matter how slow the 
quench is even for a finite system. 

The one-dimensional XXZ model is defined by the 
Hamiltonian 
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describing N spin-^ interacting via a nearest-neighbour 
Heisenberg interaction anisotropic along the z-direction, 
A being the anisotropy parameter. Here cr^, ct^ and tr^ 
are the Pauli matrices. This system is invariant un- 
der rotations around the z axis, so that the total z- 
component of the spin S^^^ is a conserved quantity. For 
time- independent couplings, the system can be exactly 
solved, by means of the Bethe ansatz [2^. If A > 1 
the system is ferromagnetic, with all spins aligned in 
the z-direction. The low-lying excitations have a gap 
l\E = 4J(A - 1) which closes for A -> 1+. For A < -1 
the system is in the antiferromagnetic Neel phase. The 
two possible ground states, differing by a traslation by 
one lattice spacing, are degenerate in the thermodynamic 
limit, i.e., for iV — > 00 their splitting is oc e"*^^, with c a 
constant. The low-lying excitations are made up by do- 
main walls separating regions with the two different Neel 
phases characterizing the ground state: they have a total 
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Figure 1: Excitation energy of the two lowest excited states 
of the XXZ model for A*' = 100 spins with open boundary 
conditions in the subspace Sfot = 0. Data obtained from 
static DMRG simulations (m = 160 with 3 target states). 
Details of the spectrum in the region close to A = 1 are 
shown in the inset. 

magnetization S^^i = 0, ±1 and a finite gap which again 
closes for A —1^. The part of the spectrum which is 
relevant for our pourposes belongs to the S^^^ — sector. 
The excitation energies of the two lowest-lying excited 
states in this subspace for a system of = 100 spins and 
open boundary conditions (OBC) are shown in Fig.[TJ In 
the whole region — 1 < A < 1 the spectrum is gapless. 
For finite sizes N and — 1 < A < 1, the gap vanishes 
linearly in ^ AE sa 2ttv/N, for all values of A (with 
different A-dependent velocities v), except for A = 1, 
where the scaling is quadratic, AE « l/N"^. Therefore 
one expects a final density of defects strongly dependent 
on whether, during the quench, the system has crossed 
this point. 

In order to study the time-dependent XXZ model we 
have to resort to numerical simulations. The results 
we present have been obtained by means of the time- 
dependent Density Matrix Renornialization Group algo- 
rithm (t-DMRG) with a second order Trotter expansion 
of H [20, [m (see [s^] for a review) . We considered chains 
up to = 200 with open boundary conditions. For the 
slowest quenches we had to restrict N to smaller values 
in order to obtain reliable results. The smallest Trotter 
time-steps were chosen to be 5t = 10"'' J, and the trun- 
cated Hilbert space dimension in the DMRG was up to 
TO = 200. We checked that in all the cases presented the 
results do not depend on the Trotter discretization 5t and 
on the DMRG-truncation of the states. 

The fully polarized ferromagnetic ground state (i.e., 
all spins up or down) is a good eigenstate for every value 
of A. It is therefore relevant to consider only quenches 
that start either from the antiferromagnetic or from the 
critical region. The anisotropy parameter A is changed 
in time according to A[t) = t/rq. In this work we con- 



sidered two different situations, i) An evolution from 
the antiferromagnetic ground state with an initial value 
of the anisotropy A^ ^ — 1 to the ferromagnetic region 
at a final value Af ^ 1; ii) an evolution from A^ = 
in the critical region to the antiferromagnetic region at 
A/ « -1. 

In order to monitor the loss of adiabaticity after the 
quench, we consider the excess final energy of the sys- 
tem relative to the ground state in the given subspace, 
conveniently rescaled. More precisely 

p . . _ {^{t)\ H{t) \^{t)) - {^Gs{t)\ H{t) \^Gs{t)) 
'°'^> {^o\H{t)\^o)-{^Gs{t)\H{t)\^Gs{t)) ' 

(2) 

where Vfo is the initial wavefunction, which we take to be 
the ground state of the initial Hamiltonian Hq = iJ(ii„), 
and '^Gs{t) is the instantaneous ground state of H{t) 
(in the subspace S^^^ — 0). The denominator normalizes 
the excess energy to the maximum possible attainable 
value, corresponding to a wavefunction \^{t)) — |^'o) 
which does not evolve at all (totally impulsive regime). 
When t — > tfjij, this quantity approaches a value -Eros = 
-£'rcs(ifin)5 which coincides with the final number of de- 
fects (apart from a constant factor) for tfin/Tg S> 1, as 
only the 2;-polarization of the spins counts, in that limit, 
in determining the final energy of the system. iJi-os nat- 
urally takes into account only the defects formed during 
the quench, and ranges from E-^cs = 1 for a totally impul- 
sive situation (the wavefunction does not evolve at all), 
to E-ccs = for a fully adiabatic evolution. 

i) Antiferro to Ferro quench - The system is ini- 
tially in its ground state at A ~ —20 (we tried differ- 
ent initial values, observing no difference), and is then 
quenched at finite rate to a final value of the anisotropy 
A = 20. The residual energy E^csit — ifin) as a func- 
tion of the quench rate 1/tq for various chain lengths 
is shown in Fig. (H^A). After a saturation region occur- 
ring for very fast quenches, -Eres obeys a power-law, with 
an exponent which is approximately 0.25 (an extrapola- 
tion to the thermodynamic limit would give an exponent 
0.251 ± 0.004, see the inset of Fig. [IJA)). The origin of 
this power-law can be undestood by means of a LZ argu- 
ment supposing that the loss of adiabaticity is en- 
tirely due to the closing of the gap at A = 1 (see Fig.[T]). 
The probability of getting excited by this point only is 
given by Pc^{tq,N) = e"''''^^^^^") , where 7 is a con- 
stant related to the slope of the two approaching eigen- 
values (ground and excited states), and AE^r = 27r^ J/iV^ 
is the finite-size smallest gap at the A = 1 point, corre- 
sponding to a spin- wave of momentum k = n/N (for 
OBC). The density of defects for large tq can be esti- 
mated by evaluating the typical length L^{tq) of a defect- 
free region (in units of the lattice spacing a), e being a 
small quantity of our choice, denoting the probability of 
getting excited. Requiring Pc-^{tq, L^) = e immediately 
implies that = TQ7'/loge~^. Consequently for the 
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Figure 2: (Color online) (A) Final excess energy after a 
quench for XXZ chains of various lengths A'^ (see legend) as 
a function of the quench rate 1/tq. The dashed lines rep- 
resent power-law fits of the data for the various A'''s. Data 
from t-DMRG simulations (m from 20 to 30, 5t from 10~^J 
to 10~^J). Inset: exponents of the power-law fits for var- 
ious chain lengths as a function of 1/A^. The dashed 
line is a linear fit to the exponent, extrapolating to 1/4 the 
thermodynamic limit. (B) Evolution of the excess energy 
of a state during quenches for A' = 100 spins in the region 
— 10 < A < 10 for various rates I/tq. Data from t-DMRG 
simulations (m — 30, At — 10"'' J). 



residual energy we get: 



Ere 



'Q 



(3) 



in very good agreement with the numerical data. Thus 
it seems that the scaling of the point A = 1 dominates 
over the rest of the critical region in a fairly wide re- 
gion of quench rates. To understand the importance of 
this point, and to see how adiabaticity is lost during the 
quench, we looked at the evolution of the residual energy, 
Eq. [21 for various quench rates, see Fig. [HB). What we 




Figure 3: (Color online) Final energy after a quench in XXZ 
chains of different lengths (see legend) as a function of the 
quench velocity 1/tq (only 'slow' quenches plotted). Data 
from t-DMRG simulations (m = 20, 5t from 10"^ J to 10"^ J). 
The dashed lines represent LZ prediction starting from the 
ground state (lower curves) or from the first excited state 
(upper curves, see text), for the various lengths. 



find is however more complicated than expected. The 
wavefunction is not frozen throughout the critical region. 
Nevertheless at A = 1 there is a clear kink in the depar- 
ture from the adiabaticity which dominates the density 
of defects (and the residual energy) in the final state. For 
the slowest quenches the dependence of the residual en- 
ergy on the quench rate crosses-over from a power-law 
to an exponential. As it has been pointed [ij, i this 
regime is described by the LZ theory. Here a good esti- 
mate of the residual energy can be obtained by assuming 
that the whole behaviour is determined by what happens 
at A = 1 where there are two excited states that become 
degenerate at A > 1. Lower and upper bounds to the 
residual energy can be simply obtained by considering 
the LZ transition probability from the ground state to 
the first excited state (lower bound) and the transition 
from the first to the second excited state (upper bound) . 
As it can be seen from Fig. [3] the actual rates lie between 
these two curves. 

ii) Critical to Antiferro quench - Another interest- 
ing situation which can be studied with the XXZ model 
is the adiabatic quench from the gapless phase, for ex- 
ample at A = 0, to a final point deep inside the Neel 
phase. In this case the critical region terminates with 
a Berezinskii-Kosterlitz-Thouless point at A = —1, and 
the applicability of the KZ theory needs to be tested. 
This situation is also relevant for the cases of strongly 
interacting bosons in one dimension driven from the su- 
perfluid to the Mott phase [3]. 

In our simulations we let A to evolve linearly from 
Ai = to A/ = —6. The residual energy Eres as a 
function of the quench rate 1 /tq is shown in Fig|3) The 
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Figure 4: (Color online) Final residual energy after a quench 
in XXZ chains of different lenght N (see legend) as a function 
of the quench rate 1/tq. The dashed lines represent a power- 
law fit of the data (apart from saturation). Data from t- 
DMRG simulations (m — 30, 5t = 10~^J). Inset: exponents 
of the power-law fits for various chain lengths A'' as a function 
of l/N. The dashed line is a linear fit of the exponents. 



results are apparently similar to the previous case. There 
is a saturation region for very fast quenches, which turns 
into a power-law for slower quenches. The exponent of 
the power-law decay is however different in this case: an 
extrapolation to the thermodynamic limit of the expo- 
nent gives in this case a value 0.78 ± 0.02 (see the inset 
of Fig. m . We believe that this is a true manifestation 
of the crossing of a critical line. Given the critical expo- 
nents of the system, a LZ (as well as KZ) treatment of 
this evolution for a single critical point would instead give 
an exponent 0.5. Differently from the previous case, all 
the gaps encountered during the evolutions have the same 
scaling and almost the same intensity, as can be seen from 
FigUI In fact in this region the system undergoes a non- 
adiabatic evolution through the critical region, which is 
not well described by the models previously proposed. 
We also note that the arguments recently put forward in 
Ref. [1^ are not applicable to the present case. 

In conclusion, the XXZ model provides a new paradigm 
to study adiabatic dynamics in many-body critical sys- 
tems. We showed in this work that depending on the 
type of quench the defect formation is dominated by the 
presence of a critical point or of a critical line. When the 
quench occurs between the antiferromagnetic and ferro- 
magnetic phase the scaling for the defect density can be 
obtained by focusing on the loss of adiabaticity at A = 1. 
In the other case we considered, when the quench starts 
from the critical region and ends in the antiferromag- 
netic region, the exponent obtained signals the crossing 
of a critical region. 
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